Location of gap nodes in the organic superconductors k-(ET) 2 Cu(NCS)2 and 
K-(ET) 2 Cu[N(CN) 2 ]Br determined by magnetocalorimetry 
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We report specific heat measurements of the organic superconductors k-(ET)2Cu(NCS)2 and k- 
(ET)2Cu[N(CN)2]Br. When the magnetic field is rotated in the highly conducting planes at low 
temperature (T ~ 0.4 K), we observe clear oscillations of specific heat which have a strong fourfold 
component. The observed strong field and temperature dependence of this fourfold component 
identifies it as originating from nodes in the superconducting energy gap which point along the 
in-plane crystal axes {d xy symmetry). 

PACS numbers: 74.70.Kn,74.25.Bt 



INTRODUCTION 



The organic charge transfer salts with general for- 
mula k-(ET) 2 X [where ET stands for bis(ethylenedithio)- 
tetrathiafulvalene] have attracted considerable interest 
because of their unconventional properties. 1 Like the 
high-T c cuprates this family of quasi-two-dimensional 
materials exhibit a low temperature superconducting 
ground state which is in close proximity to an antifcrro- 
magnetically ordered Mott insulating state. The position 
of the various members of the series in the phase diagram 
is determined by the 'chemical pressure' exerted by the 
anion X. There is considerable evidence that the super- 
conductivity is unconventional^ Power-law temperature 
dependencies observed in thermal conductivity, 2 NMR, 3 
magnetic penetration depths and recently specific heat 
measurements^ point strongly to there being nodes in 
the superconducting energy gap in certain directions of 
fc-space. In order to gain a better insight into the mech- 
anism for the superconductivity it is important to know 
the location of these nodes. 

The two most widely studied organic superconduc- 
tors are k-(ET) 2 Cu(NCS) 2 and K-(ET) 2 Cu[N(CN) 2 ]Br 
(which we abbreviate to k-NCS and K-Br respectively) as 
these have the highest superconducting transition tem- 
peratures (T c of ~ 9.5 K and ~12.5K respectively) at 
ambient pressure. Angle dependent magneto-thermal 
conductivity- measurements indicate that in k-NCS the 
gap symmetry is d x 2_ y 2 (i.e., with nodes at 45° to crys- 
tal axis). This conclusion is also supported by angle de- 
pendent tunnelling measurements,— however it is at odds 
with most theories of superconductivity such as that of 
Schmalian 8 based on spin-fluctuation mediated pairing 
(for a review see Ref. 0) which predict a d xy pairing state 
(i.e., with nodes along the crystal axes). In this case the 
nodes are in the same location, with respect to the dom- 
inant nearest neighbor antiferromagnetic coupling direc- 
tion, as in the high T c cuprates. However, Kuroki et 
alM^ calculate that a d x 2_ y 2 pairing state is often close 
in energy to the d xy and can dominate for certain model 
parameter values. 



Here we report an investigation of the location of the 
nodes in the order parameter of both k-NCS and K-Br 
using angle dependent magnetocalorimetry as a probe. 
The magnetic field was rotated in the highly conducting 
plane of the sample and oscillations of the specific heat 
were observed. The location of the maxima and min- 
ima of the fourfold component of these oscillations points 
to the order parameter having d xy symmetry. This ap- 
pears to contradict the thermal conductivity results but 
recent sophisticated theories of magneto-oscillatory spe- 
cific heat and thermal conductivity predict rather compli- 
cated phase diagrams, with the phase of the oscillations 
with respect to the nodes changing in distinct regions of 
field and temperature space. Hence, the interpretation 
of the experimental data is not straightforward and in- 
deed these two seemingly opposite conclusions are not 
necessarily incompatible. 

Angle dependent magneto-specific heat oscillations 
were first used successfully to locate the minima 
in the superconducting energy gap of YNi 2 B 2 C. 11 
Subsequently, experiments were performed in several 
heavy fermion supercondutors: Sr 2 Ru04^ CeColn 5 ^ 
Celrln 5f ii PrOs 4 Sbi2^ and UR^Si^ as well as the 
strongly anisotropic s-wave superconductor CeRu 2 ji£ In 
the simplest case, the oscillations arise from the field in- 
duced 'Doppler' shift of the energies of the quasiparticle 
states^ SE cx v s • k, where v s is the velocity of the 
screening currents (cx H) and k is the quasiparticle mo- 
mentum. Close to a node or deep gap minimum, this field 
induced shift can cause a substantial change in the pop- 
ulation of the quasiparticle energy levels. In the case of 
a simple two dimensional d-wave superconductor, if the 
field is applied along a node then only two nodes will con- 
tribute to the change in the density of states as the field 
is perpendicular to the other two. However, if the field 
is applied at 45° to this direction all four nodes will con- 
tribute, leading to an increase in the angle averaged den- 
sity of states. Hence, at sufficiently low temperature and 
field, the direction of the maxima in the density of states 
should indicate the antinodal directions.— Although sev- 
eral experimental quantities are sensitive to these oscil- 
lations, the specific heat is perhaps the most direct and 
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simplest to interpret. For example, to interpret the os- 
cillations in the thermal conductivity, scattering of the 
quasiparticles from the vortex lattice needs to be taken 
into account lit 



II. EXPERIMENTAL DETAILS 

Our samples were grown using a standard electrochem- 
ical methodic and weighed ~500/xg. The crystal orien- 
tation was determined by x-ray diffraction. Specific heat 
was measured using a purpose built calorimeter— based 
on a bare chip Cerno»^ thin film resistor. The samples 
were cooled slowly (~ 0.2K/min) to avoid any stress in- 
duced phase separation.— The Cernox was used as both 
heater and thermometer. Two different methods of mea- 
surement were used. To measure the temperature depen- 
dence of C the long-relaxation technique was used£i2i 
For the rotation studies, where much smaller tempera- 
ture excursions from the base temperature are required, a 
thermal modulation technique was usedj 22 ' 23 An ac cur- 
rent (to ~ 3-6 Hz) was passed through the Cernox and 
the signals at ui and 3w detected with lock-in amplifiers. 
The 3cj signal is inversely proportional to the specific heat 
provided that u> is selected appropriately . 22 ' 23 The ther- 
mometers were calibrated in field using a 3 He vapor pres- 
sure thermometer (below 4K) and a capacitance ther- 
mometer at higher temperature. The angle dependence 
of the thermometer magnetoresistance was measured di- 
rectly, and was found to be negligibly small. Any angle 
dependence of the addenda specific heat was checked by 
measuring a pure Ag sample at our base temperature and 
in fields up to 14 T. No change was detected to within a 
precision of 0.2 %. The thermal modulation method has 
the advantage of high resolution at the expense of a small 
(~10%) systematic error in the absolute values. 

For measurements as a function of field angle the 
calorimeter was mounted on a mechanical rotator whose 
rotation plane was parallel to the field. Sample alignment 
was done by eye and was checked using optical images. 
Using this arrangement we estimate that the field was 
kept parallel to planes to within a few (< 5) degrees. Ex- 
periments were done either by rotating the sample at low 
temperature or by heating the sample above T c after each 
rotation. In general, the two methods gave very similar 
results however there were occasional reproducible jumps 
in the data in certain field directions when the sample 
was not heated above T c after each rotation, presumably 
because of vortex pinning effects. 



III. RESULTS AND DISCUSSION 

For B ± —1AT (applied along the interlayer direction) 
both samples are in the normal state and the data follow 
C/T = 7 + /3 3 T 2 + f3 5 T 4 , with 7 = 29± lmJK^mor 1 
for K-Br and 33± 1 mJK -2 mol -1 for k-NCS in good 
agreement with previous measurements^ We do not ob- 
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FIG. 1. (Color online.) The temperature of the heat capacity 
in zero field and S x = 14 T for samples of K-Br and k-NCS 
plotted versus T 2 . The inset shows zero field specific heat data 
(with normal state contribution at B 1 ^ — 14 T subtracted.) 



serve any Schottky-like upturns even at high field, so the 
£> ± = 14T data is likely to be very close to the normal 
state C at zero fields The zero field data with the 14 T 
data subtracted shown in Fig. [T] shows the superconduct- 
ing anomalies of each sample with midpoints of 9.6 K and 
12.5 K for k-NCS and K-Br respectively. 

We now discuss the data with the field rotating in 
the basal plane. Fig. [5] shows the raw C{<j>) data for n- 
Br at several temperatures and £?H=3T which is ~10% 
of H C 2 for this field orientation 2 ^ {<f> is the angle mea- 
sured relative to the in-plane c-axis). Focusing on the 
lowest temperature data T=0.4K, we see that the to- 
tal change in C with angle is around 33% of the to- 
tal (minus addenda) or around 55% of the electronic 
term at this field (the phonon contribution at T =0.4 K 
Cphonon = 0.85 mJ mol _1 K _1 which is 60% of the total). 
The largest component has twofold symmetry, however 
to get a good fit we need also to include a fourfold term. 
We find the data are best fitted by the function 

CO) = C + \C 2 cos(0 + S 2 )\ + C 4 cos(40 + S 4 ) (1) 

which fits the data perfectly within the noise as shown by 
the residuals displayed in the figure. The fit parameters 
are C 2 = 0.54 naJmol" 1 ^ 1 , S 2 = 13 ± 2°, C 4 = 0.04 
mJinoh^- 1 , and <5 4 = 0±2°. 

The C2 two-fold term likely originates from a combi- 
nation of sample misalignment and the anisotropy of the 
Fermi surface of these compounds. If the sample is mis- 
aligned by an angle ip to the plane of rotation, then there 
will be a component of magnetic field B 1 - perpendicular 
to the planes given by / Bq = cos sin 95, where Bq is 
the applied field and <fi is the in-plane rotation angle. As 
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FIG. 2. Color online. Raw data for K-Br taken at three 
temperatures 0.4 K, 0.8 K and 1 K. The blue lines are a fit to 
a twofold angle dependence whereas the red lines include the 
additional fourfold term as described in the text. 



H C 2 is much lower for field applied perpendicular rather 
than parallel to the planes there will be a induced compo- 
nent to C(<j>) which will depend on B 1 - . For a pure d-wave 
superconductor at low field we would expect that at low 




FIG. 3. Color online. Sketch of the Fermi surface of k-NCS,— 
along with the definition of <j>- The red dots show the position 
of the nodes in for the d xy pairing state. 



fields AC{B) ~ (B , however in general the func- 
tional dependence may not be a square root (for example 
if the field is less than a scale set by the impurity scatter- 
ing bandwidth then AC ~ -B / B c2 \n(B / B c2 ) Ref. H). 
Experimentally, for these compounds for B^ > 0.1 T 
AC(B) is close to linear.— Given that this misalignment 
term must be even with respect to B, the simplest form 
is then AC(<p) = AF\ cos(0 + S)\, where A is constant 
which depends on the misalignment angle and the form 
of the function T and 5 is a constant which depends 
on the misalignment angle. For simplicity we took T 
to be a simple power law so the misalignment term is 
AC(0) = A|cos(0 + S)\ n . A least squares fit to the 
data with this in place of the C 2 term in Eq. Q] gives 
n = 1.1 ± 0.1, so we fixed n — 1 for all temperatures and 
fields. 

The crystal structures of K-Br and k-NCS are or- 
thorhombic and monoclinic respectively. The Fermi 
surface of k-NCS has been determined by both tight- 
binding 27 and first-principle 24 calculations and quantum 
oscillation measurements — and is sketched in Fig. [3] The 
quasi-two-dimensional Fermi surface would have a near 
ellipsoidal cross-section if all the ET-dimers were equiva- 
lent. However, the difference in the inter-dimer hopping 
integrals causes a gap to appear where the ellipse crosses 
the Brillouin zone boundary. The Fermi surface of K-Br is 
very similar to that of k-NCS except that there are twice 
the number of Fermi surface sheets because the unit cell 
contains twice the number of formula units due to the 
doubling of the c-axis.— 

Theoretical calculations based on spin- fluctuation pair- 
ing suggest that the pairing state has either d xy or d x i_ y i 
symmetry. For the case of d xy pairing as the gap has 
nodes along the crystal axes&2 so the angle between these 
nodes remains 90° even in the presence of the orthorhom- 
bic distortion, and maxima in C(</>) with 4- fold symme- 
try as in Eq. [1] are expected. In the semiclassical the- 
ory the field dependence of the specific heat depends on 
the Fermi velocity and gradient of the superconducting 
energy gap A near the nodes dA/ d(j>\ no d e ^ So an or- 
thorhombic distortion of the Fermi surface can produce 
a two-fold | cos</>| dependence of C(B) similar to the case 
of field misalignment. Indeed, in the field angle depen- 
dent thermal conductivity measurements of k-NCS a two 
fold term around two times larger than the four fold term 
was observed at T — 0.43 K (Ref.Q) even though the field 
was aligned to better than 0.01° . 

As the phase of the two- fold term in our measurements 
implies symmetry far from a crystal symmetry direction, 
this suggests that a sizeable fraction of it originates from 
misalignment, however the thermal conductivity results 
suggest that the intrinsic contribution is not negligible. 
In principle, in-situ adjustment of the crystal orientation 
(using for example a vector magnet as in Ref. 0) and an- 
gle detection using the magnetoresistance of the sample 
would allow us to distinguish between these two contri- 
butions however this was not possible with our current 
set-up. 
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FIG. 4. Color online. Angle dependent heat capacity data 
with the twofold terms subtracted. The angle dependent 
changes (C(4>) — C2(</>)) are divided by the angle averaged elec- 
tronic specific heat at each temperature and field C e (H,T). 
The solid lines are fits to C4 cos (40 + 84). 



In Fig. 2] we show the K-Br data with the twofold term 
subtracted and normalized to the angle averaged value of 
the electronic specific heat at the relevant temperature 
and field C e (H, T). At T=0.4K the fourfold term is now 
clearly evident. It has a relatively large peak to peak 
amplitude of ~7% of C e (H,T). Note that the in-plane 
shape of this sample is an irregular polygon and so the 
fourfold term cannot result from shape dependent vortex 
pinning effects. The phase of this term is 64 = 0±2°, so 
that the minima occur when the field is along the a and 
c crystal axes (equivalent to the b and c axes for k-NCS 
in Fig. [3]). Hence, we conclude that the nodes in the 
gap function are located along the crystal axes, i.e., the 
energy gap has d xy symmetry. 

As the temperature is increased the relative size of the 
fourfold term compared to the angle average decreases 
rapidly, and is barely discernable in the T=1K data. 
This strong decrease with increasing temperature points 
strongly to the origin of the fourfold oscillation being 
from nodal quasiparticles, rather than, for example, Vf 
anisotropy which would have a much weaker T depen- 
dence. For CcRu2, which has a strongly anisotropic 
s-wave gap (A max /A m i n ~ 5 — 10), it is found that 
C^/Ce increases with increasing T up to a maximum at 
T/T c ~ 0.16J£ Fig. 2] also shows that the magnitude of 
the fourfold term depends strongly on the applied field. 
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FIG. 5. The normalized amplitude of the fourfold component 
for both k-NCS and K-Br as a function of temperature at B = 
3T. The inset show the field dependence at fixed T — 0.4 K. 
The solid lines are fits to the nodal approximation theory (see 
text), the dashed lines in the inset are guides to the eye. 



It is maximal at i3=3T and is significantly reduced for 
5=2Tand4T. 

Data for k-NCS are also shown in Fig. [U The mag- 
nitude of the fourfold term is comparable to that found 
for K-Br, and importantly the maxima/minima are at the 
same angles, so we expect the symmetry of the gap func- 
tions to be the same. For this sample we had to introduce 
an additional term C\ cos(^> + 5\) which is comparable in 
magnitude to the C2 term to fully fit the background, 
the origin of which is unclear. The temperature and field 
dependence of the fourfold component C4 for both sam- 
ples is shown in Fig. [5] As for K-Br, Ci/C e in k-NCS 
decreases strongly with increasing temperature and de- 
creasing field. However, for k-NCS Ci/C e decreases sig- 
nificantly less as the field is increased from B=3T to 
B=4 T than in K-Br. 

In Fig. [5] we compare our experimental data to quan- 
titative predictions of the magnitude of this effect for a 
d-wave superconductor. For simplicity we use the nodal 
approximation of the Doppler shift theory which gives 
the following expression for the field/angle dependence 
of the density of states^^^ 



N(u,H,T) 
N 





f u \ 











E2 



(2) 



where the field scale Eh = 0.5ahvpj x 2 y/nfioH/^o, 
Ei = E H \ sin(7r/4 - <f))\, E 2 = E H \ cos(tt/4 - (f>)\ and 



'tf [1 + 1/(20% 



if y>i, 



(1 + 2y 2 ) arcsiny + 3j/yl — y 2 /yw, if y < 1 



In this approximation, which was shown to be in good 
agreement with more sophisticated treatments in the 
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low temperature/field limit ^ only quasiparticles in the 
nodal directions are included. This density of states was 
used to calculate the entropy S then numerically differ- 
entiated to give the heat capacity C = TdS/dT. For the 
curve in the figure we set Eh = 0.1 3\/fioH fcsT c . 

The theory and experiment show reasonable agree- 
ment considering this contains just one fitting parame- 
ter Eh- In particular, the theory explains the rapid de- 
crease in signal with increasing temperature. For re-NCS, 
using the in-plane Fermi velocity vp — 1.1 x 10 5 m/sj2i 
and penetration depth anisotropy 7a — 100, - this gives 
Eh — 0.17y7*oH ksT,, (setting the numerical constant 
a = 1). This agreement with the experimental value is 
probably fortuitous as a and the effective value of 7a de- 
pend on the details of the vortex latticed The increase 
in field between 2 and 3 T is similar to that predicted 
but the theory does not explain the observed decrease in 
higher field. However, more sophisticated treatment o 17 ' 30 
actually predict that the oscillations in C change sign for 
T > 0.ir c or for H > (0.2 - 0A)H c2 and so a strong de- 
crease of C4 with increasing field is expected near these 
transition points. 

IV. CONCLUSIONS 

In summary, the magnetic field angle dependence of 
the heat capacity of re-Br and re-NCS shows a clear four- 
fold oscillation. These oscillations are shown to arise from 
nodes in the superconducting gap situated along the crys- 
tal axis. These results imply a d xy order parameter that 
is consistent with the spin fluctuation mediated pairing 
theory of Schmalian^ and others.— At first sight these 



results seem to contradict the angular thermal conduc- 
tivity re experiments performed on re- NCS 6 which showed 
that the small fourfold component (~ 0.1% of the total 
re) was 45° out of phase with the one observed in this 
experiment (i.e., a maximum of re was observed with the 
field along the b or c axes whereas a minimum is observed 
here). A similar discrepancy between the C and re mea- 
surements was found for CeCoIns j 13 ' 32 Vorontsov et alr^ 
have shown that the oscillations in C and re both change 
sign as a function of H and T and so the CeCoIns re- 
sults could be explained because the measurements of C 
and re were conducted in different regions of the (H,T) 
phase diagram. According to the phase diagrams in Ref. 
l33l our C measurements should be well within the low 
(T,H) limit so the minima are along the nodal direc- 
tions. The re measurements were made at similar H and 
T and also should be in the low (H,T) limit with minima 
along the nodal directions, however in this case they are 
much closer to the sign switching phase boundary. Hence 
the previous thermal conductivity and the present data 
are not necessarily inconsistent and might be explained 
if there were small material dependent changes to the 
global phase diagrams of Ref. |3J. Further work will be 
required to understand this. 
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